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O ■ ABSTRACT 

(N ; 

^ , Context. The nearby spiral galaxy M81 harbors in its core a Low-Luminosity AGN (LLAGN), and appears closely related to the 

Q ' more distant and powerful AGNs seen in quasars and radio galaxies. The intrinsic size of this object is unknown due to scattering, 

' and it has shown a core-jet morphology with weak extended emission rotating with wavelength. 

I Aims. The proximity of M 81 (Z) = 3.63 Mpc) allows a detailed investigation of its nucleus to be made. The nucleus is four orders of 

. magnitude more luminous than the Galactic centre, and is therefore considered a link between Sgr A* and the more powerful nuclei 

' of radio galaxies and quasars. Our main goal was to determine the size of M 81* at a shorter wavelength thus directly testing whether 

the frequency-size dependent law oc v"" ** was still valid for wavelengths shorter than 1 cm. In addition, we also aimed to confirm 
the rotation of the source as a function of frequency. 

Methods. We observed the continuum 7 mm radio emission of M81* using the Very Long Baseline Array on Sep 13, 2002, using 
nearby calibrators to apply their interferometric observables to the target source, to increase the chances of detection. The source was 
detected on all baselines and hybrid mapping was possible. 
(~| \ Results. We present the first 7 mm VLBI image of the core of M81*, which represents the highest resolution image ever of this 

source. Modeling the interferometric visibilities with two Gaussian functions sets constraints on the angular size of its core down to 38 
I ■ microarcseconds, corresponding to a maximum (projected) linear size of 138 AU, and shows extended emission towards the NE with 

O ' a position angle of ~ 50°. A fit of one Gaussian elliptical function yields a position angle of 28 ± 8 degrees for its elongated, compact 

}~l ' structure. Combining the 7 mm size with earlier measurements at other frequencies we determine a frequency-size dependence of 

C/3 ' oc y(-OM±OM} 

Conclusions. Our VLBI imaging of M8r has clearly detected its core-jet structure, and has allowed us to estimate a size for its 
core, with a minimum size of 138 AU (~ 100 Schwarzschild radii). Our work opens the avenue for further observations of M 81* at 
I I high-angular resolution, including the monitoring of its structure, given that much higher bandwidths are currently available on the 

, interferometric networks. In particular, this would allow testing for possible proper motions of the core or of its components in the 

■ inner jet of M 81*, as well as for the speed of the detected jet components. 

. Key words. Galaxies: active — Galaxies: nuclei — Galaxies: individual: M81 — Radio continuum: galaxies — Instrumentation: 

OO ' interferometers 

; 

y—i 1. Introduction size is more than two orders of magnitude — and its luminosity 

y~~^ about four orders — greater than that of the central source in 

^ The spiral galaxy M81 (B095 1-1-633, J095532+69038, our Galaxy. Sgr A* (e.g.. Reuter & Lesch [19961) . Therefore. 

; NGC 3031, Z 333-7, UGC5318) at a distance of D =3.63 Mpc, the proximity of M 81* allows for a detailed investigation of its 

"-j (Freeman et al. 1994) harbors, together with the S eyfert 2 galaxy nucleus to study the link between our (weak) Galactic centre 

^ Centa urus A (D~3.4Mpc, e.g., FeiTarese et al. 120071 Muller and the more powerful nuclei of radio galaxies and quasars. 

H ] et al. 120111) . the closest exti-agalactic active nucleus. Hubble Indeed, X-ray and radio observations show that the ratio of the 

. . . . Space Telescope (HST) spectroscopic ob servat ions imply a (5 GHz) radio-to-X-ray (soft) luminosity for M8r can vary 

central mass of 7 x 10^ Mq (Devereux et al. 120031 ). The nucleus from Rx ~ 1.8 x 10"'' up to values of Rx ~ 3.5 x 10""^, with the 

of M 8 1 (hereafter referred to as M 8 1*) emits in radio with former value being typical of radio-quiet radio galaxies, and the 

a compact structure and exhibits both l ow-ion ization nuclear latter of radio-loud LLAGN (see Terashima & Wilson 2003), 

emission-line region (LINER; Heckman 119801) and Seyfert 1 which suggests that M8r shares properties of both kinds of 

characteristics. It thus appears closely related to the more distant objects, or may transit from one to another with time, 
and powerful active galactic nuclei (AGNs) seen in quasars and 
radio galaxies. However, M81* occurs in a spiral rather than 

in an elliptical galaxy and is relatively small and faint. In fact. As seen with VLBI (Bietenholz et al. 120001 Marti- Vidal et 

very-long-baseline interferometry (VLBI) observations (e.g., al. 120101) . the nucleus of M 81* shows a stationary feature in its 

Bartel et al. 1982) have shown that the central source is about structure — tentatively identified with the core — with a one-sided 

1000^000 AU across, depending on observing frequency, and jet to the north east. The apparent source size is ~ 0.5 mas 

that the radio luminosity of M81* is a: 10^^^ ergs"', which (~1800AU) at 8. 4GHz, and follows apower-law with frequency 

classifies M81* as a low-luminosity AGN. Nonetheless, its (0 oc v"*' **) between 2.3 GHz and 22 GHz. The sky orientation 
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of the radio structure of M 81* seems to be also frequency de- 
pendent, changing from ~ 75° at 2.3 GHz to ~ 40° at 22 GHz. 

Variability of M 8 1 ' has been reported in several occassions, 
from X-rays, both long-term (e.g., Ishisaki et al. 1 19961 ) and short- 
term (e.g., Markoff et al. i2008i) . up to millimeter (Schodel et 
al. |2007; Markoff et al. |2008] | and radio (e.g.. Ho et al. [19991 
Marti- Vidal et al. 2010) wavelengths. The turnover frequency of 
the synchrotron spectrum is in the range 150-200 GHz (Schodel 
et al. 120071 Doi et al. 120111) . In particular, VLA observations 
by Ho et al. ( 119991 ) at wavelengths of 2 cm to 20 cm in the mid 
1990s showed outbursts lasting up to three months. Marti- Vidal 
et al. (1201 Oi l has recently shown that the peaks of flux density 
are shifted with frequency when variability occurs, which can be 
interpreted as being due to opacity effects in the inner regions 
of the jet. M81* is also very peculiar, together with SgrA*, in 
its polarisation properties, since it shows significant circular po- 
larisation (Brunthaler et al. 2001), rather than being linearly po- 
larised, which is much more common among AGNs. 

VLBI observations at frequencies above 22 GHz were lack- 
ing for M81*, limiting our ability to set even tighter constraints 
on the size of M81*. Even from space- VLBI observations at 
5 GHz (Bartel & Bietenholz 120011 ). the source structure re- 
mained unresolved. Therefore, we observed this radio source at 
a frequency of 43 GHz (Al mm) using the Very Long Baseline 
Array (VLBA), with the main goal of getting the size of the 
object at 43 GHz. In turn, this would represent a direct test of 
the apparent frequency-size dependence (0 oc y^^^) observed at 
lower frequencies. 



2. Observations and data reduction 

We observed M 81* on September 13, 2002, using the complete 
VLBA (Napier [T99n i at a frequency of 43 GHz in left-hand cir- 
cular polarisation. Data were recorded with 2 bit sampling at an 
aggregate data rate of 256 Mbit s"', splitting the band in 8 inter- 
mediate frequency (IF) channels of 8 MHz each (full bandwidth 
of 64 MHz). 

The Seyfert 1 Galaxies 3C 147 and 3C 286 were used as pri- 
mary calibrators (fringe finding, delay offsets, etc.) of the ob- 
servation. We used as the phase-reference calibrator the BL Lac- 
type object B0954+658 (J0958+6533, z=0.368, 34arcmin apart 
in the sky). From UT 14:43 to UT 20:30 we cycled between 
the target source M8r (60 sec) and B0954+658 (30 sec). The 
source B095 1+699 was also used for phase-referencing test- 
ing purposes during 30min (UT 14:10 to UT 14:43) as a phase- 
reference test target. 

Without phase-referencing considerations, at the observed 
data rate, the baseline sensitivity of the VLBA for an integration 
time of ~ 60 s was of 66 mJy (1-cr) for a sub-band of 8 MHz as 
used in the observations (23 mJy for the whole band of 64 MHz). 
So, a source with a total flux of ~ lOOmJy and resolved struc- 
ture will hardly be detected by long baseline on a sub-band. 
Assuming that the source would be detected on all VLBA base- 
lines, 3-hr of integration would provide an image with a 1 cr ther- 
mal noise of 0.26 mJy beam The overall observing time was 
5.75 hr Data were correlated at the Array Operations Center of 
the National Radio Astronomy Observatory (NRAO) in Socorro, 
NM, USA. 

We followed standard procedures to a first calibration of 
the amplitudes and phases using JilVS. An a priori ampli- 
tude calibration was performed using measured system temper- 
atures and gain curves. After 2-bit-sampling digital correction, 
amplitude calibration and the removal of the parallactic angle 



phase, a single-band delay and phase offsets were calculated au- 
tomatically with the measured phase-cal values at each antenna. 
The North-Liberty VLBA station was used as reference antenna 
during the whole procedure. The Brewster antenna was flagged 
from the data due to bad performance. 

The main process of the data reduction, the search for 
group delay and phase rate calibration with the task fring in 
JilVS, had to be performed in several steps. Direct FRiNcing 
on B095 1+699, 3C 147, and 3C286 did not provide a satisfac- 
tory percentage of detections. For B0954+658 we got in a first 
attempt detections more than 80% of the time. The data for this 
source were then exported and imaged in difmap, where a few 
iterations of phase and amplitude self-calibration were applied. 
The amplitude was only corrected for an overall factor for each 
antenna. This factor was applied to all sources back in JilVS. 
The CLEAN components of the data resulting from difmap were 
also imported back to JirVS and was used as input for the task 
FRING to obtain new values for the phase and rate for the calibra- 
tor source. Those values were then transferred to the the target 
source M81*, and the JilVS task imagr was applied to these 
data to obtain a preliminary phase-referenced image of our tar- 
get source. 

Once we had a solution for the rates and delays for M81*, 
transferred from B0954+658, we narrowed the search windows 
in FRING and performed a new search for the target source, with 
a signal-to-noise satisfactory detection threshold of 3, to get a 
satisfactory number of visibilities with a high detection rate. The 
resulting data set could then be exported to difmap and hybrid 
mapping with phase and (eventually) amplitude self-calibration 
was performed on M 8 1 * . 

3. Results and Discussion 

3.1. Phase-referencing 

The phase-referenced image obtained after this process had a 
weak peak with a small offset with respect to the position 
used in the correlator. The value for the position was obtained 
by applying the JilVS task jmfit to the image obtained with 
IMAGR. The JMFIT position is therefore a - 09''55'"33?172932 
and S = 69°03'55''0610744. The uncertainty provided by jmfit 
for both values (which does not account for systematics in the 
astrometric data reduction) is 0.12milliarcseconds. This should 
not be compared directly with the phase-referenced positions to 
SN 1993 J reported by other authors, since another reference has 
been used. Different attempts to get a phase-referenced image af- 
ter removing different antennas were done, and we noticed that 
the first two hours of observations and the removal of the anten- 
nas in Hancock (system temperature, Tsys, values of above 150 K 
and relative humidity values, /i,., of ~ 40% were measured; the 
other antennas had values of ~ lOOK and dryer atmospheres) 
and Saint-Croix (T^ys values above 200 K and hr values above 
60%) from the data set had the effect of improving the phase- 
referenced image. 

3.2. Hybrid mapping 

For the data set with direct detections by fring after using 
the calibrator's solutions (that is not the phase-referenced so- 
lution and we preferred this approach due to the bad weather 
during the observations producing high tropospheric systemat- 
ics), we could perform hybrid mapping with standard proce- 
dures, involving multiple iterations of CLEAN and phase self- 
calibration. The resulting angular resolution for the observed 
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Fig. 1. Contour image of M 8 1 ' obtained by hybrid mapping on 
the self-calibrated data set. Contours are logarithmic, separated 
by a factor V2, with the lowest level set at 1.8mJy/beam, that 
is, 3 times the root-mean-square noise of the image. The peak of 
brightness is of 1 12 mJy/beam. The interferometric beam, shown 
at the bottom left, has a size of 0.361 x 0.168 mas of position 
angle -1.6°. 



(m, v)-distribution time is of approximately 0.36x0. 17 mas at a 
P.A. of~0°. 

Amplitude self-calibration was only performed once, with a 
smoothing time segment of 1 hr The amplitude corrections re- 
mained within 15% of the original calibration which also in- 
cluded the gain corrections for the calibrator B0954H-658 (see 
the previous Sect.). Applying shorter time intervals to the ampli- 
tude self-calibration did not provide satisfactory results, likely 
due to the large uncertainties in the visibilities. The result- 
ing hybrid map, produced with natural weighting, is shown in 
Fig. [U The root-mean-square noise reached in the image was 
560 yuJy beam ' . The emission to the SW is not significant and 
does not affect the model fitting results shown below. The struc- 
ture of M 8 1 * is very compact, with hints of emission towards the 
NE. This structure information is present in the closure-phases, 
and is a robust result. Indeed, any attempt to remove this emis- 
sion to the NE, while trying to get it in any other direction (by 
setting appropriate clean windows) or self-calibrating the data 
without jet emission, yielded unsuccessful results, with the resid- 
ual map always showing the need of emission towards the NE. 
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Note: The values for each bibliographical reference have been averaged in frequency. 
References: (1) Marti- Vidal et al. (2010);(2) Bartel et al. (1982); (3) Bietenholz et al. (1996); 
(4) Kellennann et al. (1976); (5) Bietenholz et al. (2000); (6) Markoff et al. (2008); 
(7) this work. 



3.3. Gaussian model fitting 

Using the calibrated data set, we model fitted the interferomet- 
ric visibilities with Gaussian functions to parametrize the source 
emission. We tried to reproduce the observed visibilities by us- 
ing four different fitting procedures: (i) a single, Gaussian, circu- 
lar function (component); (ii) a single, elliptical Gaussian; (iii) 
two circular Gaussians, and (iv) two elliptical Gaussians. The re- 
sults are shown in Table [T] (flux density, offset in right ascension 
and declination with respect to the main component, major axis 
of the Gaussian function, axis ratio, position angle of the major 
axis, and reduced likelihood parameter). Our results obtained for 
a single, elliptical Gaussian are comparable with earlier publica- 
tions (see Bartel et al. 1995). 

The fit of two circular Gaussians allowed us to quantify both 
the core and jet emission. The main, central component accounts 
for 80% of the whole emission, while the remaining ~2Q% is in 
the jet component, which has a position angle of ~ 50°. 

Table |2] shows earlier published data for M81* at different 
observing frequencies together with our own data (last row). The 
table contains the relevant parameters characterizing the source 
structure: major axis, axis ratio, position angle, and equivalent 
diameter The equivalent diameter is a measure of the area cov- 
ered by the Gaussian function. Our size determination for M 81* 
is given in the last row of the table, and corresponds to the best- 
fit to a single, elliptical Gaussian component (see Table 1). The 
uncertainties shown for the size determination at 43 GHz have 
been calculated from the response of x~ to small variations (in 
small steps of 1% for the major axis and axis ratio up to +10% 
and of 1° in a range ±10° for the position angle of the major 
axis). 

Figure |2] (top) shows the size of M 81* as a function of fre- 
quency, using two different parameters: (i) the size of the major 
axis of the source (red triangles), and (ii) the equivalent diam- 
eter of the source (blue circles). A power law fit to the major 
axis results in a dependence of size oc y-o **8±o.04 same fit to 
the equivalent diameter yields essentially the same dependence, 
size oc y-0 89±o.o8 rpjjg model predictions at 43 GHz for the two 
parameters are ~ 0.12 mas and ~ 0.082 mas, respectively, which 
are smaller than the measured quantities of 0.160 + 0.025 mas, 
and 0.087 ±0.037 mas, respectively (even with the very large un- 
certainties). If we compare the value for the size of the nucleus 
as taken from the one-component elliptical Gaussian fit (162 yuas 
in the longest dimension, with an equivalent diameter of 86 fias) 



4 



Eduardo Ros and Miguel A. Perez-Torres: VLBI imaging of M 81* at 43 GHz 




Frequency (Hz) 

Fig. 2. Size (top) and position angle (bottom) of MSI* as a 
function of frequency. We plot the data from Table 2. The two 
lines in the top panel are power-laws fit to the shown data. The 
slopes are, respectively, -0.88 + 0.04, and -0.89 + 0.08, for the 
major axis and the equivalent diameters of the source. The line 
in the bottom panel, drawn as a guide, is the function (p - (350 + 
36)° - (30 + 4)° X (log(v/Hz)). We did not use the data from 
Ref. 6 in Table 2, since this is the core component only, and 
an extended jet component, 1 mas away, was also fitted (in this 
way, is not comparable with the other values, which fit core and 
jet together). 



with the one obtained from the model fit including two circu- 
lar Gaussian components, from the brightest one (38 yuas) we see 
a margin for uncertainty in the fit. In any case, the comparison 
with the longer wavelengths is valid and we have found limits to 
the size of the emitting region in M 8 1 * at /1 7 mm. 

Therefore, the value found for the angular size of M81* is 
in the range (162-45)/ias (major and minor axis, respectively) 
down to a lower limit of 38 //as. Those values correspond to lin- 
ear sizes of (590-163) AU to 138 AU, respectively. Notice that 
for the central mass of 7 x 10^ Mq, the Schwarzschild radius cor- 
responds to 1 .4 AU. This is the most tight constraint on the size 
of M81* ever. We also note that size of the major axis follows 
very well a power-law with frequency, although our measure- 
ment may suggest a flattening of this trend, following the same 
method as earlier authors. The core brightness temperatures de- 
termined from the model fit sizes in Table |2] are in the range 
10*'" '"''' which is of the same order of magnitude as the val- 
ues for LLAGNs reported in Anderson & Ulvestad (I2005I I. and 
is therefore below the inverse Compton limit. 

If we plot the position angle of the major angle of the el- 
lipse as a function of the logarithm of frequency (Fig.|2l top), a 
clear linear trend is evident, with our data point confirming pre- 
vious results reported by Bartel, Bietenholz, & Rupen (119951) . 
among other authors. A logarithmic fit to the data points yields 
a negative trend of RA.= (350 ± 36)° - (30 ± 4)° log(v/Hz), 
implying that the source orientation rotates northwards with in- 



creasing frequency. Following this trend, the RA. of the core at 
the turnover frequency of ^ 200 GHz would be ~ 10°, this value 
being the most probable. Several scenarios can explain this ap- 
parent rotation: a strong jet bending as suggested by Marti- Vidal 
et al. (12010) model fitting; a wide angle opening in the base of 
the jet, visible at high frequencies; or even if the 'core' corre- 
sponds to radio emission of the accretion disk region. Future ob- 
servations at high frequencies with astrometric registration of the 
core region should give the answer. 

4. Summary 

We have presented the highest resolution image ever of the nu- 
cleus of M81*, and have set up a stringent constraint on the 
(projected) size of its core of 138AU, or ~100 Schwarzschild 
radii. By making use of existing size VLBI measurements for 
M81* from 1.7 up to 22.2 GHz, and adding up our measure- 
ment at 43 GHz, we find that the size of the core and jet region 
of M81* is best-fit by a frequency-size dependent power-law 
oc y-o 84±o.04^ agreement with previous results. Our 43 GHz 
data point may suggest, though, a flattening of the power-law at 
frequencies around, or above, 43 GHz. 

Our work opens an avenue for future, multi-epoch high- 
resolution VLBI observations at 43 GHz, as well as at higher 
frequencies, which will help elucidate some of the most impor- 
tant issues yet unsolved for M81*. In particular, such observa- 
tions would be very useful in constraining parameters of the ra- 
dio emission models for M 8 1 *, as well as providing valuable in- 
formation regarding the structure variability of the core-jet, e.g., 
the lifetime of the perturbations traveling down the jet. 
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